Far-infrared observations of the interstellar medium 
Editors: C. Kramer, R. Simon et al. 
EAS Publications Series, Vol. ?, 2008 



EXTRAGALACTIC H3O+: SOME CONSEQUENCES 

Floris van der Tak^, Susanne Aalto^ and Rowin Meijerink^ 



Abstract. We discuss some implications of our recent detection of ex- 
tragalactic H3O''': the location of the gas in M82, the origin of energetic 
radiation in M82, and the possible feedback effects of star formation 
on the cosmic ray flux in galaxies. 
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^ 1 Introduction 
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, Last year saw the first detection of the H3O+ molecule outside the Galaxy (jVan der Tak et al 
Using the new 16-pixel HARP imaging spectrometer on the James Clerk 
" ' Maxwell Telescope, line emission at a rest frequency of 364 GHz was detected to- 
wards the prototypical starburst galaxy M82 and the prototypical ultraluminous 
— . merger Arp 220. The derived H3O+ abundances and H3O+/H2O ratios imply 
CO ' very high ionization rates f or the dense mo lecular gas in the nuclei of these two 
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galaxies. Chemical models ([Meiierink et al. 2007) indicate that the origin of this 



high ionization is irradiation by X-rays in the case of Arp 220, whereas for M82, 
a combination of ultraviolet light and cosmic rays is needed. The high ionization 
rates of these galactic nuclei make magnetic fields more effective in retarding their 
gravitational collapse and th e formation of stars. The orig in of the irradiation in 
M82 is the evolved starburst ( Forster Schre iber et al."2003V whereas for Arp 220, 



an AGN may be needed, as claimed before by Downes & Eckart (2GG7J. Naturally, 
some questions remain, of which this paper discusses a few. 

2 Location of the H3O+ in the nucleus of M82 

The velocities and widths of the two components of the H3O+ line profile observed 
toward M82 may be used to constrain the location of the gas. Observations of CO 
emission lines toward the M82 nucleus show a flattened structure of size wSO x 20" 
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(1000x400 pc) with emission peaks on its north-eastern and south-western ends. 
The velocity field is well described by a monotonic gradient along the major axis, 
so that the struct ure is probably a rot ating ring or torus. Observations of CO 6-5 
with the JCMT ijSeaauist et al.ll2006h show a third peak in the middle, wh ich is 
even better visible in interferometric images of CO 1-0 (| Walter et al.ll2002h and 
known as the 'central hotspot'. 

The simple velocity structure of the rotating ring implies that the velocities of 
the HaO"*" components correspond to unique positions. In particular, the velocity 
of the narrow component of 14ioiio~270km s^"'^ corresponds to the inner edge of the 
north-eastern peak, which is just covered by the beam of the H3O+ observations. 
The width of the narrow component agrees with that of the CO emission from the 
north-eastern peak. 

The velocity of the broad HaO"'' component of Vhciio~220km s^^ corresponds 
to a position right between the two main CO peaks. Since the width of this 
component is much larger than the systematic velocity gradient within one JCMT 
beam, the most plausible origin of the broad component is the 'central hotspot' of 
molecular emission. This hotspot is located close to the dynamical center of M82 
and shows large line widths in several other molecular lines. The distribution of 
HsO"*" in M82 thus appears to mimic that of other dense gas tracers such as HCN 
(jMauersberger fc Hen"kelll99lh . 



3 Origins of energetic radiation in M82 

The X-ray luminosity (Lx) of M82 on arcminute scales as measu red by ROSAT 



and A SCA is ~2 x 10''^ erg/s after correction for internal absorption (jMoran fc Lehnert 



19971 ). In the central < s?ml00i,pc of M82, where HaO"*" has been detected, both 
the nuclear non-thermal component and the central thermal component may con- 
tribute to the ionization of the molecular gas; the superwind is too far away and 
its lumin osity too low for it to play a role. Observations with Chandra and XMM- 
Newton ()Zez as 2006) seem to indicate that the emission in the extended nuclear 
region is dominated by high-mass X-ray binaries (HMXBs) and so-called ultra- 
luminous X-ray sources (ULXs). The X-ray emission from a starburst is domi- 
nated by HMXBs, s uch that Lx may be used to estimate its star formation rate 



(jCrimm et al.ll2003f ). Future spatially and spectrally resolved observations of H2O 
line emission in M82 may be used to determine the H3O+/H2O ratio for both com- 
ponents of HaO^ emission and to characterize the ionization state of the molecular 
gas in M82 as a function of position in the nuclear disk. 



4 Effects of star formation on the cosmic-ray flux 

The main astrophysical interest in H3O+ is the use of the H3O+/H2O abundance 
ratio as a probe of the ionization rate of dense molecular gas. This rate is an 
important parameter for the ability to form stars, because it determines whether 
magnetic fields may be effective in supporting the cloud against gravitational col- 
lapse. The ionization rate of dense molecular clouds is dominated by ultraviolet 
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light and X-rays near their surfaces, and low-energy cosmic rays {E^l GeV) in 
their interiors. Foreground absorption often hides these types of radiation from 
direct observation, which is why H3O"'" observations are so useful. 

Recently, Pellegrini et al.l ( 2007 ) proposed a model for the photodissociation 
region M17 where the thickness of the atomic gas layer between the ionized and 
molecular components is determined by magnetic pressure rather than gas pres- 
sure. The motivation for this model is the unusually strong magnetic field mea - 
sured in M17 through the Zeeman effect on the HI 21 cm line ( Brogan et al.[l999[ ). 
Similar results w ere fo und fo r the neutral gas in front of the Orion nebula (the 
so-called veil) by I Abel et al.l (2004). The idea is that soon after the formation of 
a star cluster, starlight momentum and gas pressure from the ionized region com- 
press the surrounding neutral gas, where the magnetic pressure builds up until it 
can withstand the external pressure. See also Ferland (this volume). 

A consequence of Pellegrini's model would be that the cosmic-ray flux in the 
gas surrounding the star cluster is also enhanced, if the cosmic rays are trapped 
in the magnetic field. The efficiency of this trapping depends on the geometry of 
the magnetic field, in particular whether the field lines are open or closed. Direct 
tests of this model are not easy because the Zeeman effect is difficult to measure. 
However, an indirect test would be to look for correlations between the cosmic-ray 
ionization rates of star-forming regions with their stellar luminosity. 

Figure [1] shows recent observational estimates of the cosmic-ray ionization rate 
in dense molecular clouds plotted against their luminosity. Although the estimates 
are probably only accurate to order of magnitude, a trend does appear to be visible. 
The possible correlation cannot just be a distance effect because the ionization 
rates are local quantities derived from the ratio of two molecular abundances. 
Thus, star formation may indeed influence the local cosmic-ray flux. 

Besides variations with luminosity, the cosmic-ray ionization rate is also known 
to differ between diffuse and dense clouds. Ionization rates in diffuse gas are 
factors of 3-10 higher than those in dense clouds in the same region which are 
exposed to the same incident cosmi c-ray flux. Th i s effe ct is show n by the open 
symbo ls in Figure [U with data fr om Le Petit et al. ( 2004) for C Per, McCall et al 



(|2nn2h for Cyg 0B2, lOka et all ^00^) for Sgr A and beballe et al.l ^200(f loi 
IRAS 08572. The offset between the diffuse and the dense clouds must be due 
to pro pagation effects, and s cattering of cosmic rays off plasma waves may play 
a role ( Padoan &: Scald 2005 ). but only at low column densities {Ay ^ 10). In the 
bulk of the clouds, the lifetime r of the cosmic rays is probably limited b y energy 
losses, which scale as r 2 x 10^ {hh /300cm~^)~^ yi (jCabici et al.ll2007h . 



The authors thank Gary Ferland, Frank Israel and Sera MarkofT for useful discussions. 
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Fig. 1. Relation between cosmic-ray ionization rate and bolometric luminos- 
ity. Ionization rates are from ICaselli et af] (j2002h. IVan der Tak fc van Dishoec^ 
(|2000l). Wan der Tak et~all ()2006l) and Van d cr Tak et aD ()2008l); luminosit i es are from 
IVan der Tak et all lj2000h . lLis fc CarlstromI (jlOOi l and lSoaans fc Meiierinlj lj200'if ). The 
luminosity of the pre-stellar core L1544 is an upper limit. 
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